Universal condition for critical percolation thresholds of kagome-like lattices 
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Lattices that can be represented in a kagome-like form are shown to satisfy a universal percolation 
criticality condition, expressed as a relation between P3, the probability that all three vertices in 
the triangle connect, and Po, the probability that none connect. A linear approximation for P3(Po) 
is derived and appears to provide a rigorous upper bound for critical thresholds. A numerically 
determined relation for P3(Po) gives thresholds for the kagome, site-bond honeycomb, (3-12^) lattice, 
and "stack-of-triangle" lattices that compare favorably with numerical results. 

PACS numbers: 



Percolation is the study of long-range connectivity in 
random systems. The value of the site or bond occupa- 
tion probability where that connectivity first appears is 
percolation threshold Pc [1]. Finding exact and approx- 
imate PcS for percolating systems on various lattices is 
a long-standing problem that continues to receive much 
attention todaWe^ B i, i, H, 0, 0, 0, i, H El, H H 

All known exact pcS are for two-dimensional lattices 
that can be represented as arrays of triangular units self- 
dual in the triangle-triangle (A-A) transformation, as 
illustrated in Fig. [T] for the case of a simple triangular 
array. When this duality is satisfied, Pc is determined by 
the simple condition 0, l2l| 



p^ 



(1) 



where P3 is the probability that all three vertices of the 
triangular unit connect, Pq is the probability that none 
connect, and the prime indicates a A-A-dual system. 
The shaded triangular units can contain any collection 
of bonds, including correlated bonds which can mimic 
site percolation, connecting the three vertices. 

If, for example, the triangular unit is simply a triangle 
of three bonds, each occupied with probability then 
Pq = and P3 = p^ -\- Sp'^q, where q = 1 — and 
^ yields the bond criticality condition for the trian- 
gular lattice as q^ = p^ -\- 3p^q which has the solution 




FIG. 1: A-A duality for lattices in simple triangular ar- 
ray, (left) Shaded triangles represent any collection of inter- 
nal bonds, (right) Result of A-A transformation where blue 
(dark) triangles are the dual triangles, and form the same 
arrangement as on the left but rotated 180°. 




FIG. 2: (left) Shaded triangles in the generalized kagome con- 
figuration, (right) Result of A-A transformation, showing 
that this system is not self dual. 



Pc = 2sin7r/18 = 0.34729636 [22]. Likewise, taking a 
star of three bonds as the basic unit gives Pq = q^ -\- ?>q^p 
and P3 = p^, and ([T]) yields q^ + Zq^p — p^ or pc = 
1 — 2sin7r/18 = 0.65270365 for the honeycomb lattice 
0]. Eq. ([1]) has been applied to many other lattices 
that satisfy A-A duality, including "martini" [il. ItI. . 
bowtie [5, 2^, and "stack-of-triangle" [20] lattices, to find 
exact PcS. 

However, when A-A duality is not satisfied, then Eq. 
^ cannot be used to find Pc. For example, the A-A 
transformation for the kagome lattice is shown in Fig. [21 
and it can be seen that, while the lattice can be broken 
up into non- touching shaded triangular units, the A-A 
transformation gives a different lattice altogether, and 
so the self-duality condition is not satisfied. Likewise, 
site percolation on the honeycomb lattice, which can be 
represented as bond percolation on the kagome lattice 
with all three bonds correlated (see Fig. |3(a)[ ), is also 
non-self-dual. 

Nevertheless, for any system that can be broken up 
into identical disjoint isotropic triangular units, Pc must 
be determined by a unique condition that depends only 
upon the connections probabilities Pq and P3 of the tri- 
angular units. In this paper we consider lattices of the 
kagome form, as shown in[4l[d), and investigate the corre- 
sponding relation between P3 and Pq. The kagome form 
includes several unsolved lattices of interest as discussed 
below. While we can't find exact thresholds for these 
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FIG. 3: (a) Site percolation on the honeycomb lattice (red) is 
equivalent to bond percolation on the kagome lattice (shaded 
triangles) with all three bonds occupied, or all vacant, (b) Ba- 
sic unit for analyzing site-bond percolation on the honeycomb 
lattice in the generalized kagome framework. 



lattices (indeed, they are likely insolvable), we can make 
very precise predictions on their values and unify their 
study. 

First we consider the "double honeycomb" lattice, 
shown in Fig. St^b), which is of the kagome form and is 
the one exactly soluble lattice of this form. It can be con- 
structed by replacing each bond of a honeycomb lattice 
(Fig. [ll^a)) by two bonds in series, which implies that its 
Pc is the square root of the Pc for the honeycomb lattice: 



p 



-2sin7r/18 = 0.80790076 



(2) 



For this lattice, which we indicate by a star, we have 



p 



Sq^'^p'' = 0.09652861 (3) 
: 0.52731977 (4) 



where = 1 — p^. Note, Eq. ([T]) is far from being satis- 
fied. 

Next, generalizing the considerations in [4], we develop 
an approximate linear relation between P3 and Pq for all 
lattices of the kagome form, that is exact at the point 
(Pq, P3*). Consider the systems shown in Fig. [H In (c) 
we replace all the up-stars of (b) with general shaded 
triangular units with a given net connectivity Pq and 
P3. This produces a generalized "martini" configuration, 
which falls under the general triangular class of Fig. [H 
with connectivities (as follows from the diagram in (c)): 



M) 
^3 



Po + 3P2(^*' + 2gV) + Psiq""^ + V) 
PsP*' (5) 



Eq. ([T]) then yields the exact criticality condition for sys- 
tem (c): 



P3 = Ps^b{Po-P^) 



(6) 



where b = 1/(2 -p*) = 0.83885634. As a final step, we 
hypothesize that Eq. ([6]) represents an approximation to 
Pc of the "full" kagome system with both up and down 
triangles shown in Fig. H^d). The justification is that in 
going from (a) to (b), we replaced one set of stars by 
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FIG. 4: Steps in the derivation of the linear relation Eq. 
(|6|): (a) the honeycomb lattice, (b) double- honeycomb forms a 
kagome class of lattice, (c) all up-stars replaced by triangular 
units, forming martini configuration satisfying A — A duality, 
(d) remaining stars replaced by triangular units, forming the 
kagome configuration. 
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FIG. 5: Plot of Ps - [Ps* + b{Po - Pq*)] vs. Pq - Pq using 
data of Table HIl showing deviations from Eq. 0. Points 
are numerical data, and the curve is a plot of Eq. ©. The 
locations of some specific systems are also shown. 



shaded triangles satisfying (|6]), and the system remained 
at criticality. Now we replace the second identical set of 
stars by the same shaded triangles, and we expect that 
the system remains close to criticality. 

In Table [J we compare the predictions of the linear re- 
lation (|6]) with the numerical results for several systems. 
The (linear) estimates are found by putting the cor- 
responding expressions for Pq and P3 into Eq. (|6|) and 
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TABLE I: Results of pc and Po, P2 and P3 for various system. ^Determined by Eq. 
[23, ""Ref. [9], ^Ref. [14], ^this work. Pq, P2 and P3 are calculated using pc (cubic) 



©, determined by Eq. {|9)), ^Ref. [4], ^Ref. 



system 


Pc (linear)" 


Pc (cubic) ^ 


Pc (numerical) 


Po 


P2 


P3 




U.OU 1 i7Uvj 1 u 






U.Uc/UOZ/OU-L 




'^97'^1Q77 


(3, 122) 


0.74042118^ 


0.74042081 


0.74042195(80)^ 


0.10045606 


0.12297685 


0.53061341 


kagome 


0.52440877^'^ 


0.52440516 


0.52440499(2)^ 


0.10757501 


0.11861544 


0.53657867 


honeycomb (site) 


0.69891402 


0.69702981 


0.69704024(4/ 


0.30297019 





0.69702981 


00 subnet 






0.628961(2)^ 
0.625365(3)^ 


0.09652861 


0.12538387 


0.52731977 


subnet 4 


0.62536437 


0.62536431 


0.09823481 


0.12433811 


0.52875085 


subnet 3 


0.61933204 


0.61933180 


0.6193296(10)^ 
0.6008624(10)^ 


0.10016607 


0.12315455 


0.53037028 


subnet 2 


0.60086322 


0.60086202 


0.10402522 


0.12078995 


0.53360494 



solving numerically for p. For the kagome lattice, we use 
Po=q\ Ps=p^^3p\. (7) 



<7 , n=P 

For the (3, 12^)-lattice (shown for example in Ref. [4]) we 
use 



Po = 1 - 3p^ - 3p^ + 6//^ + 3/ - 4//^ 
Ps = 3//2 - 2//2 . 



(8) 



For site percolation on the honeycomb lattice, Pc = P3 = 
1 - Po, and Eq. ^ yields explicitly Pc = ^/[p'^'^i^ - 
p'^)] = 0.69891402. The agreement between (linear) 
and numerical results is especially good for systems where 
Po is near Pq . 

To test the behavior of P3(Po) over a more complete 
range of values, we carried out new simulations using the 
gradient percolation method [25|, [2^ on a general kagome 
systems. We fixed Po = 0,0.5,0.1,0.15, and 0.25 and 
allowed P3 to vary linearly in the vertical direction, with 
the estimate of the critical value found as the fraction 
of Pa-triangles in the frontier. We considered systems 
of different gradients and extrapolated the estimates to 
infinity to find the values of P3 given in Table HIl 

In Fig. Owe plot the difference between the measured 
P3 and the predictions of Eq. (|6|) as a function of Po for 
these systems. The first derivative at Po = Pq appears 
to be zero, which would imply that Eq. (|6]) represents 
the exact linear term in the behavior of P3 vs. Po — ^0 • 
The numerical data also suggests that (|6]) gives an upper 
bound for P3(Po) for all Po. Fitting the data to a cubic 
equation, assuming that P3(Po ) = b exactly, we find 

Ps = Ps + b{Po - Po*) + c(Po - Po*)' + d{Po - Po*)' (9) 

with c = -0.05987 and d = -0.1038. This curve fits 
all the data points P3 within ±10~^. The results of 
using this equation to predict Pc are shown in Table [J 
under the heading "cubic", and all are within the ex- 
pected error of about ±10~^, and more accurate as Po 
approaches Pq . For the kagome case, our prediction 
Pc = 0.52440516 compares favorably to the recent pre- 
cise result 0.52440499(2) of Ref. [H (which appeared af- 
ter our analysis was complete) and the previous value 
0.5244053(3) [27]. 




FIG. 6: Lattices with subnets 2, 3 and 4 (left to right) 



We next apply our general relation for P3 vs. Po to 
get very accurate PcS for a class of lattices in which each 
triangle of the kagome arrangement contains a "stack-of- 
triangles" as shown in Fig. [6l In Ref. [io] the similar 
stack-of-triangles were studied in a regular triangular ar- 
rangement, and explicit expressions for Po and P3 were 
found by exact enumeration for these three subnets. We 
can use those same expressions to analyze the subnets on 
the kagome lattice as well. For subnet 2, we have [lO] 



Po 

P3 



q-' + y 



■P 
(10) 



with q = 1 — p. For subnets 3 and 4, see Ref. [20[. 

We insert these expressions for Po and P3 into Eqs. 
(j6|) and (|9]) to find the linear and cubic estimates for 
Pc. The resulting values are shown in Table H along with 
results of numerical simulations. For subnets 3 and 4, the 
predictions of (|6|) and especially ([9]) are expected to be 
very accurate, because Po is so close to Pq , and indeed 
the precision of the numerical simulations is not high 
enough to see the difference between these predictions 
and the actual values. 

As seen in Table IH the quantities Po , P2 and P3 ev- 
idently approach the double-honeycomb values Pq , P2 
and P| as the mesh of the subnet gets finer. This is be- 
cause the triangular units in the fine-mesh limit can be 
effectively represented by a star of three bonds, with the 
central site in this star representing the supercritical "in- 
finite cluster" in the central region of the triangular units 
(iot . The set of these stars creates the double- honeycomb 
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TABLE II: Results of simulations for P3 and P2 = (1 - Po - 
P3)/3 for general kagome systems as a function of Po; val- 
ues are accurate to about 10~^. These data are plotted in 
Fig. O Also shown are the equivalent site-bond probabilities 
Ps and pb calculated from Eqs. (|12|) . The third row is the 
double- honeycomb system and the final row represents site 
percolation on the honeycomb lattice [14]. 



Po 


P2 


Ps 


Pb 


Ps 





0.1846972 


0.4459084 






0.05 


0.1539432 


0.4881704 






0.0965286 


0.1253839 


0.5273198 


0.6527036 


1 


0.1 


0.1232560 


0.5302320 


0.6583497 


0.9926153 


0.15 


0.0926739 


0.5719784 


0.7405771 


0.8974788 


0.2 


0.0622208 


0.6133375 


0.8242773 


0.8195766 


0.25 


0.0319205 


0.6542385 


0.9091230 


0.7547482 


0.3029598 





0.6970402 


1 


0.6970402 



expression we found, Eq. (|6]), appears to be exact to first 
order, and the simulation results shown in Fig. [5] sug- 
gest that that expression provides upper bounds to Pc 
for these systems. We conjecture that this is indeed the 
case. The numerically refined cubic relation of Eq. [9] al- 
lows very accurate thresholds to be predicted for a wide 
variety of systems, and an explicit expression for the crit- 
icality condition of site-bond percolation on the honey- 
comb lattice to be written. 

This work was supported in part by the U. S. National 
Science Foundation Grant No. DMS-0553487. 



lattice, so the Pi are the same as the double-honeycomb 
values. Furthermore, the probability Poo, corner of con- 
necting from a corner to the central infinite cluster at crit- 
icality must be identical to the double-honeycomb bond 
threshold, p'^. Thus, we can find Pc for the infinite net by 
running simulations of growing clusters from the corner 
of a single large triangular system, and adjusting p until 
Poo,corner(p) = P^ This yiclds Pc{oo) = 0.628961(2). 

Finally, we note that a realization of the general 
kagome system for Pq > Pq* is given by site-bond per- 
colation on the honeycomb lattice, as represented in Fig. 
3(a) For the site-bond basic unit of Fig. |3(b){ we have 



Po 
Ps 



1-Ps +p.[(l - v^)' + 3(1 - v^)'V^] 



3/2 



which can be inverted to yield: 



Pb 



3P3 



2P. 



(11) 



Ps = Ps/P^ (12) 



In Table [m we list the values of and Ps that correspond 
to the measured values of P3(Po). We can also put Eq. 
(pT]) into Eq. (|6|) and simplify using Eqs. (|3|) and (|4j) to 
find an approximate expression for the critical line on the 
Ps-Pb plane: 



(13) 



Pb{l-B{^-p-)) 



where B = p'^ / — p'^'^) . We can improve upon this 
relation by using the cubic function of P3(Po) given in 
Eq. ([9]); this adds the additional terms C {^/ph — p'^)'^ + 
D{^Jph— p^Y above formula, where C = 9p'^^(2 — 

p*)V(3 -p*^)^c = -0.0460682 and D = -0.01681. 

In conclusion, we have shown how the notion of a 
unique relation between P3 and Pq, first studied in the 
context of self-dual systems 0, [2l| , extends to the non- 
self-dual kagome configuration. The approximate linear 
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